Integrable semi-discretization of the coupled 
nonlinear Schrodinger equations 



Takayuki Tsuchidafj], Hideaki Ujino and Miki Wadati 

Department of Physics, Graduate School of Science, 
University of Tokyo, 

Hongo 7-3-1, Bunkyo-ku, Tokyo 113-0033, Japan 
Received 29 April 1998, in final form 25 November 1998 

Abstract. A system of semi-discrete coupled nonlinear Schrodinger equations is 
studied. To show the complete integrability of the model with multiple components, we 
extend the discrete version of the inverse scattering method for the single-component 
discrete nonlinear Schrodinger equation proposed by Ablowitz and Ladik. By means of 
the extension, the initial-value problem of the model is solved. Further, the integrals of 
motion and the soliton solutions are constructed within the framework of the extension 
of the inverse scattering method. 
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1. Introduction 

There has been a surge of interest in the family of nonlinear Schrddinger (NLS) equations 
because of its many applications to various kinds of physical phenomena. In the 
remarkable papers [[I], 0] Zakharov and Shabat solved the NLS model, 

■dq d 2 q 

Ot OX 1 m j\ 

dr d 2 r 

1 »T ~ "sT~2 + 2r <? r = 
ot ox* 

with r = =pg* by means of the inverse scattering method (ISM). After another success of 
the ISM for the modified KdV equation ||, Ablowitz, Kaup, Newell and Segur ||] unified 
a class of soliton equations by employing various time dependences of the scattering 
problem. Their formulation is called the AKNS formulation. 

A number of authors have studied extensions of the AKNS formulation and 
presented many models, which are integrable by the ISM [^|, 0. Among such models, a 
system of coupled nonlinear Schrddinger (CNLS) equations 

Mj d 2 gj ? ^ _ n 



dt dx 2 ,. i 

}2„ 



dr~ d rj 



J = l,2,...,m (1.2) 



k=l 



is particularly remarkable in describing diverse physical phenomena [?| || |[| [T^, [U], 



I|, 0, g|, 0, |T|, 0, 0, 0. Manakov § considered the two-component CNLS 
equations ( (|1.2|) with m — 2) with = — g* (j = 1,2) as a model for propagation 
of two polarized electromagnetic waves and applied the ISM to the model for the first 
time. Interest has recently focused on the two-component CNLS equations in studying 
explicit solutions [H], [11], the stability of solitary waves [|14| and interactions between 
solitons in birefringent optical fibers [T^, [B| from a physical point of view. 

Very recently, the authors proposed a new extension of the ISM and solved the 
coupled modified KdV (cmKdV) equations j20| 



~0f + 6 ( £ C ^u k ) ^ + ^ = z = 0, 1, . . . , M - 1 (1.3) 



in the self-focusing case, which had been investigated by alternative approaches |21| , |22 



By a transformation of variables, this model is cast into a new coupled version of 



the Hirota equation p3| , which describes wave propagation in optical fibers, including 
higher-order effects. 

On the other hand, some discrete versions of the ISM have been constructed and 



applied to some discrete models E4L Eq, Q% Gq, 29 1. Among those models, the 



semi-discrete nonlinear Schrddinger (sd-NLS) equation found by Ablowitz and Ladik 
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28| 



+ g n _i - 2q n ) - g„r n (g n+1 + g n _ x ) = 



.dq n . ( 

i^r - (r n+ i + r n _! - 2r n ) + r n g r , 



;i-4) 



n+l 



+ r 



n-1 J 







has been studied extensively because of its simplicity and physical significance. They 
solved ( |1.4|) under the rapidly decreasing boundary conditions, q n ,r n — > at n — > ±oo. 



The model (1.4) was also solved under other integrable boundary conditions [30, 31. 32 

The sd-NLS equation has attracted researchers to studies on various subjects, such 
as nonlinear lattices in condensed matter physics |£J, phase plane patterns J34|, breather 
solutions [35], Backlund transformations []36[| , numerical experiments and homoclinic 
structure [37], [38], , the dynamics of a discrete curve [[5J and surfaces , Hamiltonian 
structure and the classical r-matrix representation j|2[ [43], |4| and the quantization of 
the model pE2"f . 

In an analogous way to the continuous theory, it is natural to consider a 
generalization of the sd-NLS equation ( |1.4| ) with multiple components, namely, 



dq { J ] 

i 4^ + (q ( n 



Ot 

■ n 



0) j_ Aj) 
+1 "r q n -i 



2tf ) )-£? 



(*)-.(*) 

n n 



(J) I Jj) > 

+1 "r 5n-l> 



r Cf) 
n+l 



r 0') 
'n-1 



fe=l 
m 



2r«) + £ 



(*)„(*) . (r U) 
n Hn v n+l 



+ r, 



(3) > 
n-1/ 











1,2, 



m. 



1.5) 



k=l 



We call this model the semi-discrete coupled NLS (sd-CNLS) equations. The model 
is expected to be important in various applications, e.g. numerical simulations of the 
CNLS equations ( |1.2| ). Hisakado [|5| showed that the system ( |1.5| ) is connected with 
the two-dimensional Toda lattice. An iV-soliton solution was obtained by Ohta [|4"6]| . It 
is noted that another scheme of integrable semi-discretization of the CNLS equations 



was reported by Merola, Ragnisco and Gui-Zhang [[47 



In [48| , the authors proposed a new extension of the discrete version of the ISM by 
Ablowitz and Ladik. Applying the extension, they solved the initial-value problem of 
the semi-discrete coupled modified KdV (sd-cmKdV) equations, or the coupled modified 
Volterra equations |45| , 

duV M - 1 



dt 



1 



J2 C JkU 
j,k=0 



(?)„(*) 
n "n 



It 



(<) 
n+l 



"n-lj 



0,1, 



M-l 



1.6) 



under some appropriate conditions. A systematic procedure for constructing 
conservation laws and multi-soliton solutions was also given [ESI. Related results are 
obtained by means of Hirota's method [46 , 49| . 

In the present paper, we use a transformation of variables 



v (2j-2) + iv (2j-l) 



c ' n 



. v (2i-2) + b (2i-l) 
^n 1 ^n 



1,2, 



;i-7) 



which cast the sd-CNLS equations (|1.5|) into the sd-cmKdV equations ( |1.6|) with 



a- 



it 



(0 



and M = 2m. We pull back the transformation (|1 . 7| ) to the 
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level of the Lax representation and give an explicit Lax pair for the sd-CNLS equations 
( |1.5| ) for the first time. Following the method of |[48| , we can solve the initial- value 
problem of the sd-CNLS equations ( |1.5| ) with = —q%'* under the rapidly decreasing 
boundary conditions, q^p — > at n — > ±00. Explicit forms of conserved quantities and 
the iV-soliton solution are also given within the framework of the ISM. 

The paper consists of the following. In section 2, we introduce a Lax pair for the 
semi-discrete matrix NLS equation. Considering a reduction to the sd-CNLS equations, 
we obtain the Lax formulation and conservation laws for the sd-CNLS equations. In 
section 3, we perform the ISM for the sd-CNLS equations with directing our attention to 
the transformation ( |1 . 7|) . The initial- value problem is solved and the iV-soliton solution 
is given. The last section is devoted to discussions. 

The main idea of the paper is based on a matrix representation and some properties 
of the Clifford algebra, whose elements are anti-commutative. The proof of relations 



used in the paper is given in [4S 



2. Lax representation and conservation laws 

2.1. Lax pair for the semi- discrete matrix NLS equation 

We begin with a set of auxiliary linear equations 

* n+ i = L n m n V nit = M n m n . (2.1) 

Here \l/ n is a 2/-component column vector, and L n , M n are 21 x 21 matrices. The 
compatibility condition of ( ^.1| ) is given by 

L nit + L n M n - M n+1 L n = O. (2.2) 

We call L n and M n the Lax pair and ( |2.2j ) (a semi-discrete version of) the zero- curvature 
condition, or simply, the Lax equation. Let us introduce the following form for the Lax 
pair: 





' 






+ 









o 

F 2 R n 



F\Q n 

O 



o o 

O F 2 



zF 2 



(2.3) 



' u 
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+ z 



O lQ n 

iF 2 R n „ 1 F 1 O 



+ 



-iQ n F 2 R„-iF 1 + iHx 
O 



O 

iR n FiQ n -iF 2 + \H 2 



1 

+ - 
z 



O -iFi<5„_ 1 F 2 
-iR n O 



1 

+ ^7 



o o 

O -il 

zQn — \F\Q n -\F 2 

-^1 + R n F 1 Q n ^iF 2 + H 2 

where z is the spectral parameter which is time-independent. / is the / x / unit matrix, 
Q n and R n are / x / matrices. The constant matrices Fx, F 2 , Hi and H 2 are assumed to 
be Hermitian and satisfy 



z 2 I — Q n F 2 R n ^iFx + Hx 
zF 2 R n _xFx — ^R n 



(2.4) 



(Fx) 2 = {F 2 f = I [Fx,Hx] = [F 2 ,H 2 



O. 



(2-5) 
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into (12.2), we obtain 



Here [ • , • ] denotes the commutator. Substituting Q2.3| ) and 
a set of matrix equations 

iQn,t + Fi(Qn+i + Qn-\)F2 + H\Q n — Q n H 2 — F\Q n +\FiR n Qn — QnFLnF\Qn-\Fi = O 
iR n ,t — F 2 (R n+ i + Rn-i)F\ — R n Hi + H 2 R n + F 2 R n+ iFiQ n R n + R n Q n F2R n -iF\ = O. 

We call this model the semi-discrete (sd-) matrix NLS equation. The integrable model 
(2.6) with Fi = F 2 = I, Hi = —I,H 2 = I was found by Ablowitz, Ohta and Trubatch 



ii 



2.2. Conservation laws 

In this subsection, we present a method to construct local conservation laws for the sd- 
matrix NLS equation ( |2.6|) , which is a discrete version of the method in the continuous 
theory p0 |. We start from an explicit expression of (|2.1|) , 







Fin S n 




' *ln " 




' ' 
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B n 






. *2n+l . 




Tn F 2n 




. $ 2n . 
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C n 






. *2n . 



(2.7) 

where all the entries in vectors and matrices are assumed to be / x / square matrices. 
Introducing an I x I square matrix T n by 

r n EE ^2n^Tn ( 2 -8) 

we can show the following relations from (|2.2|) and 



Fj 



In 



A n +i — (S n T n + Fi n )A n (S n T n + Fin) 1 



+ B n+ iT n +i — (S n T n + Fin)B n T n (S n Tn + Fi n ) 

Tn+l = (T n + F 2n T n )(Fi n + s n T n ) ■ 
Taking the trace on both sides of (|2.9| ), we obtain 

tr{log(S n r n + Fm)} t = tr(A n+ i + B n+1 T n+ i) - tr(A r 
Assuming the form of L n as ( |2.3| ), we have 



B n T n ). 



Fin — zi 7 ! F 2n — —F 2 S n — FiQ n T n — F 2 R n . 



Then ( |2.11| ) and ( |2.10| ) are transformed into 

{tr log(/ + -Q n T n )) t = tr(A»+i + B n+1 T n+ i) - tv(A n + B n T r 



1 (2. 
(2.10) 

(2.11) 
(2.12) 

(2.13) 



"Z-QrXn = QnF 2 R n -lFi H QnF^n^L^Qn-lTn^ijFi — (Q n ^n)Fi (Qn-l^n-l) Fl . (2.14) 

Z 

Note that Q2.13j ) has the form of the local conservation law. This suggests that 



tr{log(J + QnTn/z)} is a generator of the conserved densities for ( |2.6| ). We substitute 
the expansion of Q n T n with respect to 1/z, 



OO 

QnX~ n = ^2 z 2j-l fl 



(J) 



(2.15) 



3=1 
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into ( 2.14 ). Then we obtain a recursion formula for fw\ 



/W = Q n F 2 R n -iFi5j tl + Q n F 2 Q-\ft 1 1) F 1 - £ f^F 1 ft 1 k) F 1 j = 1, 2, . . . . (2.16) 

k=l 

Equation Q2.16Q yields f%\ for instance, 

f [ n ] = QnF 2 R n -lF 1 

fn = QnRn-2 — QnF 2 R n -lQ n—1 

F 2 R 

n—2 ■ 

We substitute ( |2.15| ) into tr{log(J + Q n T^ /z)} and expand it with respect to 1/z, 

tr{log(/ + !/<') + 1/f + + ...)} = tr{l/« + - i(tf>)'] +...}. 

Thus, the first two conserved densities given by this expansion are 
= tr{/«} = tviQnF^^} 

= tr|Q n i?„_ 2 — Q n F 2 R n -iQ n -iF 2 R n ^ 2 — -(Q n F 2 R n -iFi) 2 ^. 



Similarly, we expand Q n T n with respect to z, 

CO 

QrJ^n = Z J &i ' 

Substitution of ( |2.18| ) into ( |2.14| ) yields a recursion formula for gftf\ 



(2.17a) 

(2.176) 
(2.18) 



(fc) 



,0'-*) 



j = 1,2,... .(2.19) 



fc=i 



From formula (|2.19|) , the first three of the coefficients g%' are given by 

QnRn 



9 P 



g^ = -Q n F 2 R n+l F x {I - Q n R n ) 

9n^ = — QnRn+2{I ~ QnRn) + {QnF^n+lF^ 2 (I — Q n R n ) 
+ QnRn+2FlQ n+ lR n+ \Fi 

\I QnRn) ■ 

We substitute (|2.18| ) into tr{log(J + Q n T^ /z)} and expand it with respect to z, 
tr{log(/ + ^ + z 2 # + z 4 ^) + •■•)} 

= tr{io g (/ + g W) + + g^r 1 + * 4 [g n 3) (i + gi l) )- 1 - + g^r 1 } 2 } +•••}■ 

Thus, the first three conserved densities in this expansion are 

J<°> = tr{log(J + = tr{log(J - Q n R n )} (2.20a) 

4 1} = ^{g [ n\l + ^i 1 ^ 1 } = H-QnFA+iFj (2.206) 

4 2) = tr[ 9 S\i + g^r 1 - \{g n 2 \i + gWY 



trj— Q n Rn+2 + + „ {.QnF 2 R n+ \Fi) 2 } . 



(2.20c) 



The generator of the conserved densities, tr{log(/+<5 n rn/ z )}, is shown to be related 
with a time-independent subset of scattering data defined later (see the appendix). 
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2.3. Reduction of the Lax pair and the conservation laws for the semi-discrete coupled 
NLS equations 

In this subsection, we show a reduction of the sd-matrix NLS equations to the sd-CNLS 
equations. We recursively define 2 m ~ 1 x 2 m ~ 1 matrices F{\ F^, H[ m \ H2 , 
and by 



(i) 



1 Fr 



(1) 



1 H\ 



(i) 



■1 H. 



(i) 



(m+l) 



(m) 



-F 



(m) 



(m+l) 



(m) 



rr(m) 1 r 



-,(m) 



(m) 



iZ 



(m+l) 



iff - /■ 



2 m-l 



2m-l 



0? = ^ ^=r« 

g(m+l) = 



„(m+l) j P(m) 



R{ rn+1) 



R (m) g("»+ 1 )/ 2m - 1 

r^ +1 )j 2m -i -Q^ 



(2.21) 
(2.22) 

(2.23) 

(2.24) 

(2.25) 
.(2.26) 



Here J 2 m-i is the 2 m 1 x 2 m 1 unit matrix. It is readily seen that (|2.5|) is satisfied. For 
the matrices defined by (|2.21|) -( 2.26 ), we can prove the following relations: 



(2.27) 



fl| m) QW - Qi m) ^ m) = -2F 1 (m) Q^F 2 (m) 

-R^H[ m) + ^ m) i?i m ) = 2F 2 (m) J R( l m )F 1 (m) 

by induction. Then substituting R^\ etc. into Q n , R n , etc. in the sd-matrix 

NLS equation Q2.6|) , we obtain the sd-CNLS equations 



at 1 n+1 



fc=i 











j = l,2,...,m. (2.28) 



For instance, the Lax matrix L n for the two-component sd-CNLS equations ( (|2.28| ) with 
m = 2) is given by 

~ z g« g( 2 ) - 
-z _rW r (D 

1/z 



r (l) fl (2) 
' n Hn 



-q ( n ] 1/2 



(2.29) 
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In what follows, we set r«' = —q^* (j = 1,2, ... ,m) and consider a self- focusing 
case of the sd-CNLS equations, 



■ n . . (gB 1 + 9£i-2g? ) ) + Ek? ) l 2 '(ffSi + rf-i) = ° J = l,2,...,m. (2.30) 
In this case, the relation (|2.27| ) for Q^ 1 ' and R^" 1 ' becomes 

m 

Q^itf = R^Q^ = - E • (2.31) 

In addition, a simple relation between Q^ m ^ and holds, 

J RM = -QMt (2.32) 



where the symbol f denotes the Hermitian conjugate. The relations ( |2.31|) and (|2.32| ) 
play an essential role in the ISM. 

Because the sd-CNLS equations are given as a reduction of the sd-matrix NLS 
equation ( |2.6| ), the results in section |2.2| assure the existence of an infinite number of 
conservation laws for the sd-CNLS equations. Explicit forms of the first four conserved 
densities for the sd-CNLS equations ( |2.28j ) are given by 

log(l - 5>S°rW) (2.33) 
i 

e*(-l) B (ff£Ui fc) - q^lUi) Vj,* 
e- 4B (-l) B (rSJi»"? ) " r«rj»i) V j, k 

&ri fc) + « 1 Vj,k (2 - 34) 

/ j Hn+X'n l^Hn 'n+X 

j 3 



(J) r (3) 



\ l / , 1n+X'n+XI / Aln+2' n ^ In 'n+2) ^ \Z-j\"n+X ' n In ' n+X) J / , Hn+X' n+X / , 1; 
j 3 3 3 j 

(2.35) 

\ L Z^y«+2' n+2) \ L / , Hn+X' n+X) / Aln+V n In ' n+3) 
3 3 3 

_ f-i _ >p„(i) J3) \ . \^(J3) (j) _L_ n U) r (j) \ . I\^( n (j) r (j) _ n (j)j3) \ i X^fJi*) r V) _ Jj) r (j) U 
\ L / j Hn+X' n+X) /_^\ L ln+2' n * In ' n+2) \/ AUn+X' n In 'n+X) T / Aln+2' n+X Hn+X ' n+2) j 

3 3 3 3 

l^ 1 / , ¥71+1 ~n+l I \ / A1n+X r n In ' n+X) / , Hn+2 r n+2 T / ,Wn+2'n+l yn+l'n+2j 'n / 

3 3 3 3 3 

+ "{E^rf - q^l)} 3 + £(&^ _ . £ gjfl^ . £ ^) r 0). (2 . 36 ) 

j j j j 

A straightforward calculation shows that all the entries in ( |2.34j ) are conserved densities. 
It is noted that the lower two densities in (2.34) correspond to the conserved densities 



qjTk and J2j Qj r j f° r the continuous CNLS equations (|1.2| ). 
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It is also remarkable that higher than second conserved densities for ( 2.28|) have 
expressions such as 

i-l 



fc=l j j 

which split into two independent densities with simpler structures, 

i-l 

(1 — Qn+k r n+k) ' Qn+i r n + ' ' 



(i>2) 



(z>2) 



and 



fc=i 



i-l 



H( 1 _ Qn+kr n +k) ■ q 



nTn+i 



+ ■ 



k=l 



(i>2) 

in the single-component case (cf equation ( |1.4| )). This fact implies that a direct recursion 
formula of the conserved densities for the sd-CNLS equations might be, if it exists, very 
complicated. We have obtained the first four conserved densities concisely by use of the 
recursion formula of the conserved densities for the sd- matrix NLS equation ( |2.6| ). 

3. Inverse scattering method 

In this section we investigate the scattering and inverse scattering problems associated 
with the 21 x 21 (I = 2 m ~ l ) matrix (JT 



zFi FtQ n 
F 2 R n ~Fl 



*2n 



(3-1) 



to solve the sd-CNLS equations (|2.30|) . Here and hereafter the superscripts (m) of F\ 



(m) 



Fr\ and are often omitted for convenience. To simplify the analysis, we 

consider a gauge transformation 



e i f n (F 1 



\n -\H\t 



Then the scattering problem ( |3.1|) is changed into 



" $ln+l " 




" zl 


Qn 






. $ 2n+l . 




_ Rn 


ij 

2 




. $ 2« . 



where the transformed spectral parameter and potentials are 
z = ze^ Q n = i n (F 1 ) n e' iHit Q n e iH2t (F 2 ) n R n = (-i) n (F 2 ) n e-' lH *R n ^\F, 
The constraints Q2.32| ) and (|2.31|) lead to those for Q n and R n , 

Rn Qn QnRn RnQn ~ 



-El^l 2 -^-^/. 

We assume the rapidly decreasing boundary conditions, 

Q n ,R n -^0 as n — > ±oo. 

Considering the time dependence of the scattering data, we can solve the initial-value 
problem of the sd-CNLS equations ( [2.30 ). Some of the main ideas in the following are 



(3.2) 

(3.3) 
(3.4) 
(3.5) 
(3.6) 



an extension and a modification of the analyses for the matrix KdV equation [p2| , the 



matrix mKdV equations and the cmKdV equations |2(| (see |48| for details) 
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3.1. Scattering problem 

Let &^\z) and 3>^(z) be matrix solutions of ( |3.3| ) composed of 21 (= 2 m ) rows and 
I (= 2 m_1 ) columns. We introduce the following matrix function of and & 2 ^: 

^ n [$( 1 ),$( 2 )] = $i 1 )t(i)$( 2 )(^). (3.7) 

This satisfies a recursion relation 

= (/ - QMW^,^] = p n W n [$M, $ (2) ] (3.8) 

where p n is defined by 

m 

(3.9) 



Pr, 



i+^ = i+Ei^' ) r- 

We introduce Jost functions 4> n , <ft n and ip n , ip n which satisfy the boundary 
conditions 



and 



I 

o 

o 
I 



o 
-I 

I 

o 



as n 



-oo 



as n — > +oo. 



(3.10) 



(3.11) 



Here O and I are, respectively, the I x / zero matrix and the I x / unit matrix. We can 
show that (fi n z~ n , ip n z n are analytic outside the unit circle (\z\ > 1) on the complex z 
plane, and that 4> n z n , ip n z~ n are analytic inside the unit circle (\z\ < 1) on the z plane, 
when Q n and R n go to O sufficiently rapidly as n — > ±oo. We assume the following 
summation representation of the Jost functions ipn and ijj n : 



Vv= E z- n 'K{n,n') i> n = E z n 'K{n,n') 



(3.12) 



where X(n, n') and i^(n, n') are ^-independent column vectors which consist of two I x / 
square matrices, 



K(n, n') 



Ki(n, n') 
K 2 (n,n') 



K(n,n') 



Ki(n, n') 
K 2 (n,n') 



We substitute ( |3.12| ) into ( |3.3| ). Equating the terms with the same power of z, we obtain 



and 



Ki(n,n) = O 

oo oo 

K 2 (n, n) = - RiQiY 1 = U p^ ■ I 

i=n i=n 

Q n K 2 (n, n) = -Ki (n, n + 1) 
K 2 (n,n) = 

oo oo 

Ki(n,n) = - Q l R i y l = ]J p, 1 ■ I 

i=n i=n 

RnKxin, n) = -K 2 (n, n + 1). 



(3.13a) 
(3.13b) 

(3.13c) 

(3.14a) 
(3.14b) 

(3.14c) 
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Since a pair of the Jost functions <p n and n , or ip n and ip n , forms a fundamental system 
of the solutions of the scattering problem (|3T3"|), we can set 

<j> n {z) = Mz)A(z) + Mz)B(z) (3.15a) 
Mz) = Mz)B(z) - i(^(4 (3.15b) 

Here the coefficients {A(z), A(z), B(z), B(z)} are n-independent I x I matrices which 
are called scattering data. 

To derive the formula of the ISM rigorously and concisely, we assume that Q n 
and R n are on compact support. The result is, however, valid for larger classes of the 
potentials Q n and R n . Using the asymptotic behaviors of the Jost functions (|3.10|) , 
( p. 11| ) and the relation (|3.8|) , we obtain 



A{z) = WJtf^] A(z) = -WoofoM] (3.16) 

and 

1 1 1 - - 1 00 

A\-)A{z) + B\ -)B{z) = A\-)A{z) + B\-)B{z)= R p n ■ I (3.17a) 

Z Z Z Z n=-oo 

A^(-JB(z) = B^-JA(z). (3.176) 

z z 

The expressions (|3.16| ) show that A(z) and A(z) are, respectively, analytic outside the 
unit circle (\z\ > 1) and inside the unit circle (\z\ < 1). 

3.2. Gel'fand-Levitan-Marchenko equations 



Multiplying A(z) 1 and A(z) 1 from the right to ( |3.15a| ) and ( |3.15b| ), respectively, we 
have 

(t) n {z)A{zY l = $ n {z) + tP n (z)B(z)A(z)- 1 (3.18a) 
Mz)A(z)- 1 = -ip n (z) + ^ n {z)B{z)A{z)-\ (3.186) 
We substitute (|3.12|) into the right-hand side of fl3.18c| ) and operate on both sides 

— I dz z- m - 1 (m > n) 

where C denotes a contour along the unit circle \z\ = 1. It should be noticed that 4> n z~ n 
and A(z) are analytic outside the unit circle C, \z\ > 1. The inverse of A(z), i.e. A(z)^ 1 , 
is given by 

where A denotes the cofactor matrix of A. We assume that 1/ det A(z) is regular on the 
unit circle C and has 2N isolated simple poles {z±, z%, . . . , Z2n} in \z\ > 1 (see ( |3.27| ) for 
the reason why we choose the number of poles to be 2N). We set 

Joo in = lim <p n z~ n A(z)~ l 
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and use the residue theorem. After some computation, we arrive at the discrete version 
of the Gel'fand-Levitan-Marchenko equation, 



(n'+m)— l 



(3.19) 



(3.20) 



K(n, m) + K ( n > n')F{n' + m) = Joo,nS n , m (m > n). 

n'=n 

Here F(n' + m) is defined by 

1 r 2N 
F{n' + m) = — d> B{z)A{z)- 1 z-^ n ' +m) - l Az + V C jZ 
2ixi jc j =1 

where Cj is the residue matrix of B(z)A(z)~ 1 at z — Zj. 
Similarly, we operate 

— / dz z m ~ x (m > n) 
2m Jc v ~ ' 

on both sides of (|3.18fl) substituting (|3.12|) . As has been mentioned previously, (p n z n 

and A(z) are analytic inside the unit circle C, \z\ < 1. We assume that 1/ detA(z) is 

regular on the unit circle C and has 2N isolated simple poles {z%, z 2 , ■ ■ ■ , z 2 n} in \z\ < 1. 

We set 

J ,„ = KmJ) n z n A{z)~ l 

z— >0 

and use the residue theorem. Finally, we obtain the counterpart of the discrete Gel'fand- 
Levitan-Marchenko equation, 



K(n,m) — ^2 K{n,n')F{n' + m) = —Jo,n$n,m ( m > n )- 



Here F(n' + m) is defined by 



1 



2N 



Fin' + m) = — i B(z)A(z)- 1 z n ' +m - 1 dz - V C k z^ +m -\ 
2vri Jc ^ 

The matrix is the residue matrix of B(z)A(z)~ 1 at z — z k . 
From ( |3.13b|) and Q3.14b|) , it is natural to set 

00 00 
K(n, m) = K,(n, m) TT (I — RiQi)~ l = n(n, m) p^ 1 (m > n) 



(3.21) 



(3.22) 



oo 



i=n 
00 



K(n,m) = K(n,m)Y[(I - QiRi) 1 



K(n, m 



\[pi l (m>n). 



Here «;(n, m) and /«(n, m) are column vectors whose elements are / x I square matrices, 



Ki{n, m) 
K 2 (n,m) 

In particular, «(n, n) and R,(n,n) are given by 



R 2 (n,m) 



«(n, n) 



O 



«(n, n) 



O 
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Due to ( |3.13| ) and ( |3.14| ), the potentials Q n and R n are given by 

-/d(n,n + l) = Q n (3.23) 

- R 2 (n,n + 1) = R n . (3.24) 

In terms of k and R, the Gel'fand-Levitan-Marchenko equations ( 3.19|) and ( [3.21|) for 
m > n are rewritten as 
O 



K(n, m) + 



F(n + rn) + ^ n(n,n')F(n' + m) 

n'=n+l 



o 
o 



K[n.m\ 



I 

O 



o 
o 



[m > n) 



[m > n). 



(3.25) 
(3.26) 



F{n + rn) — ^ n(n,n')F(n' + m) - 

n'=n+l 

It should be noted that the scattering problem (|3.3| ) gives the symmetry properties 
of the scattering data. For instance, we have 

det A(z) = det A(-z) det A(z) = det A{-z) (3.27) 

which means that the eigenvalues Zj, z^ should appear as 'positive- negative' pairs. 
Further, we have 

B(z)A(z)- 1 = -B(-z)A(-z)- 1 B(z)A(z)- 1 = -B(-z)A(-z)- 1 . 
Therefore, we can simplify the forms of F and F as 

f 2F R (n + rn) rn = n + 2j - 1 
b [n + m) = < ' 
v J \ O m = n + 2] 

1 r N 
Fr(u) = —J c B(z)A(z)- 1 z- n - l dz + Y.C J zj n ~ 1 



(3.28) 



J>1 



and 



F(n + rn) 



2Fn(n + rn) m = n + 2j — 1 



F R (n) = — B{z)A(z)- L z 

Z7T1 JCh 



O m = n + 2j 

N 

] -"- 1 dz- y £C k 2%- 1 . 



J>1 



fc=l 



Here (7r denotes the right-half portion of the unit circle contour C. 

The symmetry properties of F and F give rise to those of k and R. From ( |3.25| ) 
and (|3.26| ), we obtain 

{Km(n, m) m = n + 2j — 1 
O m = n + 2j 

{i<>2R(n, m) m = n + 2j — 1 
O m = n + 2j 

Considering the above symmetry properties, we obtain the simplified Gel'fand-Levitan- 
Marchenko equations for k 1r and R 2 r, 

Km(n, m) (m > n, m — n = odd) 



K 2 {n,m) 



J>1 
J>1- 



(3.29a) 
(3.296) 



2F R (n + m)-A ^ ^ K 1R {n, n")F R {n" + n')F R (n' + m) 



(3.30) 



n'=n+2 n"=n+l 
n'— n=even n"— n=odd 
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K2r{u, m) (m > n, m — n = odd) 

oo oo 

= -2F R {n + m)-A £ £ R 2R {n, n")F R {n" + n')F R {n' + m). (3.31) 

n'=n+2 n"=n+l 
n'— n=even n" — n=odd 



3. 3. Time dependence of the scattering data 

Under the rapidly decreasing boundary conditions (|3.6| ), the asymptotic form of the 
Lax matrix M n for the sd-CNLS equations ( |2.30|) after the gauge transformation ( |3.2| ) 
is given by 

z 2 I O 



M n = g n M n g n 1 + g nt g 1 



O 4/ 



as n — > ±oo. 



We define time-dependent Jost functions by 



J n — ~rn 



e ~ 



O 



z n e z2t = <f) n e-' - 



O 



as n — > — oo. 



From the relations 
we obtain 



1 



(3.32) 



4> n ,t = ( M n ~ Z I)<p n 4>n,t = (M n ~I)<Pn- 

Z A 

We put the definitions of the time-dependent scattering data, 
M z ) = Mz)A{z, t) + i) n {z)B{z, t) 

4> n (z) = Tp n (z)B(z,t) - 1p n (z)A(z,t) 

into ( |3.32| ). Then taking the limit n — > +oo, we obtain the time dependences of A, 
BA' 1 , Cj and A, BA^ 1 , C k . They are given by, respectively, 



A(z,t) = A(z,0) 



B(z,t)A(z,t)- 1 = B(z,0)A(z,0)- 1 e 
and 

A{z,t) = A{z,0) 



C i (t) = C i (0)e 



(3.34) 



B(z,t)A(z,t)- 1 = B(z J 0)A(z,0)- 1 e^ z2 -^ 



H C k (t) = C k (0)e K ^ . 



(3.35) 



The above results give explicitly time-dependent forms of F R (n,t) and F R (n,t) for odd 
n. 
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Thanks to the constraints R n = —Qn and Q n R n = RnQn = 
additional relations besides (13.27] ) and (|3.28|) . The first additional relation is 

detA(z) = {detA(\)r 

z* 



a n I, we have some 



(3.36) 



which is proved in This relation restricts the numbers and the positions of the 

poles of 1/ detA(z) and 1/ detA(z), i.e. 



N = N z k = -. 

4 



Due to ( |3.17£|) , we have the second additional relations, 
B(z)A(z)- 1 = {B(z)A(z)- i y (on |*| = 1) 



Ck 



Cu ■ 



: 2 k 



(3.37) 

(3.38a) 
(3.38b) 



The relations (|3.37|) and ( [3.38D give a relation between F R (n,t) and F R (n,t), 

F R {n,t) = F R {n,t)\ (3.39) 

In order to make the ISM applicable to the sd-CNLS equations, we have to take account 
of the internal symmetries of the potentials Q n and R n . According to ( |2.25| ), ( |2.26| ) and 
the potentials Q n and R n are defined recursively by 



Qi l) = iV^) 



(m+l) 



m = (-ir e - 2it rP 



Q(m) i n e 2it q n m+1) I 2 m-i 
;_i)n e - 2i M m+1 )/ 2 m-i -R^ 



R 



(m+l) 



_\\n -2\t (m+l) j 

ii e r n 'i 2 m-i 



i e q n 'j 2 m-i 



If we set 



v (2i"2) + i (2j-l) 
n 1 n 



-i) n e" 2i V^ 



J = l,2, 



(3.40) 



2m-l 



(3.41) 



_ v (2i-2) + iv (2j-l) 

and for m > 2 are written as 

2m- 1 

Qt ] = 4 0) H + E 4 fc) e fe ^ = -vS>l + E «?W 
fc=i fc=i 

Here II is the 2 m_1 x 2 m ~ 1 unit matrix, which has been denoted by I%m-i. Because of 
= —q$'*, v£> and v n k > should be real. Substitution of (|3.41|) into the relations (|3.5| ) 
yields the following important relations for 2 m_1 x 2 m ~ 1 matrices {e{\: 

{ ei ,e,} + = -25 M I (3.42) 



o fc - -e fc . (3.43) 

Here {•,•}+ denotes the anti-commutator. It is stressed that expressions ( |3.41| ) with 
conditions ( |3.42j ) and (|3.43|) are only a sufficient (i.e. not necessary) condition of ( |3.5D . 
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Considering the above symmetries of the potentials and for m > 2, we 
can show the following properties of the scattering data. 
Proposition 3. 4. 1 

(1) The determinant of A(z) and the determinant of A(z*) are related by 

det A{z) = {det A{z*)}* . (3.44) 

Thus the poles of 1/ det A(z) in \z\ > 1 appear as pairs situated symmetric with respect 
to the real axis. Therefore, we replace 2N in section 3. 2 with 4iV and choose the values 
of 2N poles in \z\ > 1, Rez > as 

, - - Wi J = l,2,...,iV (3.45) 

where dj > 1, < 6j < ir/2 for 0j 7^ 0. The conditions (|3.45|) should be interpreted as 
follows. If 9j = 0, the corresponding pole does not need its counterpart. The values of 
the remaining 2N poles in \z\ > 1, Rez < are given by 

Z2N+k = —Zk k — 1,2, ... , 2N. 

(2) The reflection coefficient B(z)A(z)~ 1 is expressed as 

2m- 1 

B{z)A(z)- 1 = r®l+ E r(k)e k- (3-46) 

k=l 

Here r^ ^ and r( fc ) are complex functions of z and t which satisfy 

r (°)( z *) = r (0) (z)* r (fe) (z*) = r (fc) (z)*. (3.47) 

(3) The residue matrices {C\, C 2 , ■ ■ ■ , C 2 n} are expressed as 

2m- 1 

C ! y_ 1 = C fl+ E <f 

J = l,2,...,iV (3.48) 

^ = cf *1 + E cf % 
fe=i 

with complex functions of time £, c^- and cf' . 

The statements (l)-(3) are proved essentially in the same way as in H3| and therefore 
their proofs are omitted. 

Taking account of the above conditions, we obtain explicit expressions of Fn(n,t) 
and F R (n,t) for odd n, 

1 r 2N 
F R (n, t) = — B(z)A{z)- 1 z- n - 1 dz + E C^" 1 

1 „ 2m— 1 

— / {(r (0) 2 _n - 1 +r (0) *z n " 1 )I+ V (r^^-™" 1 +r (fc) ^ n ~ 1 )e fc )dz 



27ri Vc^ 

JV 

+ E{( c f ^ n_1 + cf E (^f zj 11 ' 1 + # ) *^" T *" 1 )c*} (3.49) 



fc=i 

JV 2m- 1 



3=1 k=l 
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F R (n,t) = F R (n,ty 

i „ 2m— 1 

= _L / (( r (0)^-n-l +r (0)^n-l )I _ £ ( r (*)«-»-l+ r (fc)*^-l) efc }dz 
27T1 VCt/ji L ^ J 

AT 2m-l 

+ E{( c f 2 7 n_1 + cf %*- n - 1 )l - £ (cf zp-i + c f — i) efc } (3.50) 

j=l k=l 

where Cjj R denotes a contour along the quadrant (upper-right portion) of the unit circle 
C. We see that the coefficients of 3 and {e^} in ( |3.49| ) and ( |3.50| ) are real. Thus, a 
pair of F R and —Fr is expressed in the same form as ( |3.41| ), as is expected from the 
viewpoint of successive approximations for the Gel'fand-Levitan-Marchenko equations. 

Because B(z)A(z)~ 1 and Cj depend on t as ( |3.34| ), the time dependences of r^°\ 
r^ and cf , are given by 



r^(z,t) = r {0 \z,0)e^ z2 -^ t A k \z,t) = r^(z, 0)e _ ^-^)* (3.51) 



cf\t) = cf (0)e ' 3 ' ^ cf (t) = cf (0)e~ V J . (3.52) 
Combining ( |3.30| ) and ( |3.31| ) with ( |3.39| ), we arrive at 
K 1R (n,m;t) 



oo oo 



= 2F R (n + m,ty K ^( n ^ n + 2 h~ 1; ^(271 + 2/ 2 + 2Zi - 1, t)F fl (n + 2l x + m, t) f 

z 1= i i 2 =l 

(3.53) 

R 2 R(n,m;t) 

oo oo 

= - 2F fi (n + m, t) - 4 £ £ ^(n, n + 2/ 2 - 1; t)F*(2n + 2/ 2 + 2l x - 1, t) + F fl (n + 2/ x + m, t) 

h=i h=i 

(3.54) 

for m > n, m — n = odd, where F R (n, t) is given by ( |3.49| ). 

Now the initial- value problem of the sd-CNLS equations ( p. 30 ) can be solved in the 
following steps. 

(1) For given potentials at t = 0, Q n (0) and R n (0) which are expressed as (|3.41|) , we 
solve the scattering problem (|3.3|) , and obtain the scattering data {-E>(z)v4(z) _1 , Zj, Cj} 
or, more concretely {r^°\ r^ h \ Zj, cf , cf }. 

(2) The time dependence of the scattering data is given by (|3.34j ) (or (|3.51|) and ( |3.52j )). 

(3) We substitute the time dependence of the scattering data into the Gel'fand-Levitan- 
Marchenko equations ( |3.53|) and (|3.54j) . Solving the equations, we reconstruct the time- 
dependent potentials, 

Q n (t) = -K 1R (n,n + l;t) 
R n {t) = -R 2R {n,n+ l;t). 
This solution proves directly the complete integrability of the sd-CNLS equations ( |2.30| ). 
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To construct soliton solutions of the sd-CNLS equations, we assume the reflection-free 
condition, i.e. B(z) = B(z) = O on \z\ = 1. Then, F R (n,t) and F R (n,t) for odd n are 
given by 



2N 



F R (n,t) = E^'W^"" 1 = Cj(0)e 

i=i 



2N 



F R (n,t) = -J2C k (t)z%- 1 C k (t)=C k (0)e 



k=l 



To solve ( p.30| ) (or equation (|3.53|) ) with (|3.55|) and Q3.56Q , we set 



2N 



K lR {n, m; t) = E PkCkit)^- 1 (m-n = odd). 



fc=i 



Substituting ( p.57| ) into ( |3.30| ) or ( |3.53| ), we have 



27V 2N = 



^2 



1=1 j=l Z j ( Z j Z k)( Z j Z l) 

In terms of a matrix S whose elements are defined by 



- S rP l Ci{t)C j {t)=-2I. 



(3.55) 
(3.56) 

(3.57) 
(3.58) 



2N = 



^2 



Sik = hkl — 4E(~ ) —} .■■>.. ■> :~rr 

P^Zj) (z] - z 2 k )(z] - zf) 



2N 



1 



^zfzf\z)z^ -\){z)zf -\) 
equation ( |3.58| ) is expressed by 

I Sn ••• S12N \ 

( P 1 P 2 ■ ■ ■ P2N 



Ci{t)C 3 {t) 

Ci{t) ] C j{t) l<l,k<2N 



-2(J_L-__L)- 



(3.59) 



2N 



\ S2NI ■ ■ ■ S2N2N J 

Similarly, we solve (|3.31| ) (or equation ( 3.54Q ) with ( |3.55|) and Q3.56j ). Substitution of 



2/V 



K 2R {n, m; t) = E PjCj{t)i 



-(n+m)— 1 



3=1 



into (|3.31| ) or (3.54 ) gives 



2N 2N -2n+2 i 

P.-4VV zh 

3 ~2n / 7 2_;2U 7 2_z2\ 

1=1 k=l I \ Z l z k)\ z j z k) 



(m — n = odd) 



PiQMCkit) = -2L 



(3.60) 



(3.61) 



Using a matrix S, 



2N ^2n+2 



S lk = 8 ltk I-4,J2- 



W + 4E 



£ zf" (z?-z])(zi-z]) 
1 

2n^*2n/^2^*2 



Ci{t)C 3 \t) 



^zf\z}zf - \){z\zf - \) 



-d^C^tY l<l,k<2N 



-2(1 1 ■■■ I ). 

2N 



(3.62) 
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we rewrite ( |3.61| ) as 

( S\\ ■ ■ ■ S\2N \ 
( Pi P 2 • • ■ P2N ) : ' • . : 

V S2N1 ■ ■ ■ S2N2N J 

Equations ( |3.59| ) and ( |3.62j ) are readily solved. Thus the iV-soliton solution of the sd- 
CNLS equations ( |2.30| ) is given by 

QM(t) = i n {F 1 ) n e- iHlt Q^e iH2t {F 2 ) n 
= - K 1R (n,n + l;t) 




1 

1 1 

C2(ty- 



\ 



2n+2 



* 2n+2 



Z- 



(3.63a) 



2 A' 



RT\t) 



-i) n (F 2 ) n e- iH2t R^e iHlt (F i y 
R 2R (n,n + l;i) 




*1. 



c 2 (t) 



~2n+2 
z 2 



c 2N {t)- 



2n+2 



Z, 



(3.63b) 



2iV 



Strictly speaking, equation (|3.63|) includes breathers besides solitons. In order to 
extract pure soliton solutions, we assume that each soliton seen in J2j \ln (*)| 2 nas a 
time-independent shape. By calculating an asymptotic behavior of the tails of solitons 
at n — > +00, we obtain the corresponding necessary conditions 



r** r** r** — r** — 

( - J 2j-\ ( - J 2j — ^2j^2j-\ — < ^2j < -^2j-1 — ^2.3-1^2 



o 



j = l,2,...,iV 



'2j-V^2j — ^2j^2j-l — ^2j^2j-l — ^2j-1^2j 

on the residue matrices. Equation ( p.64|) is translated explicitly into 

2m-l 

./ 1.2 V. 



(3.64) 



i=0 



(3.65) 



As an example, we write down a pure one-soliton solution. Choose N — 1 and set 



-w+ie 



W > 



-1 1 2m— 1 

1 -t 1 M°)*ir ( fc )* 

z l k=l 



2n 



2m- 1 



k=l 
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C 2 =cf ] *l 



3^0 



E_(K) * 
ci e fc 

fc=i 
sinh 2W 



\ 



2m- 1 

2 E I^Wp 



(3.66) 



Then, from (|3.63|) we obtain 



Q { ™\t) = sech{2nW/ + 2(sinh2W/cos2#)t + o }- 



sinh 2W 



2m-l 



x 2 E I^WI 2 

■ {(7 1 (o)e 2i ^ ne+ ^ cosh2Wsin2e * ) *-'' + (7 2 (o)e~ 2i ^ ne+ ^ cosh2Wsin26 '* ) *-''} (3 67a) 

Rt\t)= -Q™(tf. (3.67b) 

It is straightforward to show that Q3.67| ) can be expressed as ( 3.41|) with real coefficients 
v$(t) of I and {e^.}. Thus we have checked in terms of the inverse problem that the 
conditions ( |3.41[ ) or consequently ( |3 .5| ) are satisfied under proposition 3. 4. 1 and the 
conditions ( |3.64| ), in the case of the one-soliton solution. 
By introducing a new set of constants by 

a _C^-2) (0)+i ^-l) (0) 

^^W-icf^O) 

we obtain a simplified expression of the pure one-soliton solution of the sd-CNLS 
equations ( |2.30| ), namely, 
(jw.s , r „ „, „/ . , „„, , , sinh2W 



1,2,. ..,m 



q W(t) = sech{2nW + 2 (sinh 2 W cos2#)t + O }- 

E(hf + I&i 2 ) 

3=1 

^ a e 2i{n(0- f )+(cosh 2W sin 20-l)t} _|_ ^* e -2i{n(0+f )+(cosh 2W sin 20+l)t} j 

z = 1, 2, . . . , m. 



\ 



(3.68) 



Here the condition ( |3.65|) for iV = 1 is cast into the orthogonality condition, 

m 

E OiPi = o. 



1=1 



The soliton solution ( ft-68 ) exhibits a novel property as a solution of NLS-type equations. 
Because there are two carrier waves in one envelope soliton, the shape of soliton observed 
in \q^\t)\ 2 periodically oscillates in time. It is observed for ( |3.68|) that the summation 
of |<Z„ (t)| 2 with respect to components, i(= 1,2, ... ,m) 

m 

\q%\t)\ 2 = sinh 2 2W sech 2 {2nW + 2 (sinh 2W cos 28)t + O } 

i=i 

has a time-independent shape, as is expected. This fact suggests that the conditions 
( [3.641 ) are necessary and sufficient to give pure soliton solutions even in the iV-soliton 
case. 
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The continuum limit which reduces ( |3.68| ) into the one-soliton solution of the 
continuous CNLS equations can be seen as follows. We denote by h the (dimensionless) 
lattice spacing of the sd-CNLS model. We rescale t by 



1 



and set 



hui(x, t) 

71 

~ 4 " 

A = o. 



x 







W = hi] 

OH = h^i 

If we take the continuum limit h 
transformed into 

Ui(x, t) = 2r] sech{2?7£ + 8rj^t + 0o} 



nh 



0, the one-soliton solution (|3.68| ) with ( |3.66| ) is 



-2i£z-4i(£ 2 -»7 2 )t 



\ 



1,2, 



m 



El7i 

3=1 



with 



2 V 



\ 



This is indeed the one-soliton solution of the continuous CNLS equations, 



, dm d 2 Ui 



dt 



+ 



dx 2 



+ 2E 



1,2, 



m. 



(3.69) 



,771) IS 



It is noteworthy that either a, = {i = 1,2, ... ,m) or /3j = (i = 1,2, 
necessary for us to take the continuum limit. This reflects the fact that the pure iV- 
soliton solution of the sd-CNLS equations (|2.30|) includes more arbitrary constants than 
the iV-soliton solution of the continuous CNLS equations (|3.69|) . It has its origin in the 
oscillation of solitons seen in each component \q$(t)\ 2 . In this sense, the structure of 
the pure soliton solution may be more similar to that of the cmKdV equations (|1.3j ) 
rather than to that of the CNLS equations ( |3.69| ). 

Ohta j|6| obtained an iV-soliton solution for the sd-CNLS equations in the Pfaffian 
representations. The iV-soliton solution ( |3.63j ) with the constraints ( |3.64 ) contains more 
parameters than Ohta's iV-soliton solution. Therefore, it is reasonable to conjecture that 
our solution reduces to Ohta's by a particular choice of those parameters. 

We wish to show that ( |3.63a| ) and ( |3.63b| ) are expressed as 

2m-l 

Q i n\t)=v^{t)t+ V vl k \t)e k 



2^ V n 

k=l 
2m- 1 



^l m) (*) = -^ 0) (*)l+ E «Sf°(*)e* 



(3.70a) 
(3.70b) 



fe=i 



without using the products of {e,i} such as e^ej, CiCjCk- Further, we expect to 
be real in fl3.70a| ) and ( 3.70b|) . As mentioned previously, '( [3.70D — > (|3.5|) ' holds. Thus, 
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equation ( |3.5|) is automatically satisfied if we can show ( |3.70| ). So far we have proved 



either ( |3.70a| ) or ( |3.70b| ) only for m = 2. However, the result for the one-soliton solution 



implies that both of ( P-70|) are simultaneously satisfied under proposition 3.4. 1. 
4. Discussions 

We have investigated the semi-discrete coupled nonlinear Schrddinger (sd-CNLS) 
equations. The analysis from the ISM point of view is given for the first time in this 
paper. A previous paper |48] dealt with a new extension of the semi-discrete version of 



the ISM proposed by Ablowitz and Ladik to solve the semi-discrete coupled modified 
KdV (sd-cmKdV) equations (|1.6|) . In this paper, we have developed the extension 
with the help of the transformation ( |1 . 7|) to solve the model considered with arbitrarily 
multiple components. 

We should comment on the Lax formulations and the ISM for both the sd-cmKdV 
equations and the sd-CNLS equations, which may spotlight mysterious structure of 
discrete soliton equations. 

First, there is an essential difference between the form of the L n -matrix for the sd- 
CNLS equations and that for the sd-cmKdV equations in the case of m > 2 (M > 4). 
It seems that the scattering problem for the sd-CNLS equations associated with the L n - 
matrix, e.g. ( |2.29| ), does not agree with the scattering problem for the CNLS equations 
|8|. |53[ in the continuum limit (z = e~~^ h+1,p , h — > 0). The situation is not observed 



for the continuum limit of the sd-cmKdV equations. This difference can be understood 
by considering the continuum limit of the semi-discrete coupled Hirota equations. A 
detailed explanation will be reported in a subsequent paper. 

Secondly, it is obvious that (|2.6f) does not generally allow us to assume the reduction 
R n = ±Qn, because of the order of the products. To consider the reduction R n = ±Ql, 
we should impose the additional restriction Q n R n = R n Qn = scalar on Q n and R n . 
Thus, both ( j2.31|) and ( p. 32 ) play a crucial role in our theory, which is peculiar to the 



discrete theory with multiple components. 

Thirdly, as for the consistency of the iV-soliton solution, it is difficult to prove both 
of equation ( 3.70|) from the restrictions on the scattering data ( |3.45|) -( j3.48|) . In the 



continuous theory ( |3.70| ) is easily proved at least for 2 x 2 matrices and R( 2 >. 



However, even in this case, it is not so easy to prove both of ( |3.70|) in the discrete theory. 



Fourthly, due to the relation ( 3.44Q , we always have to consider pairs of poles 



as scattering data in the ISM. The constraint gives the novel structure of solutions, 
'two carrier waves in one envelope soliton'. This situation is closely related to the 
high internal symmetries of and R^ 1 ' which may give a rich variety of dynamical 
behaviors in the multi-field soliton systems. As an influence of this fact, we need to 
assume additional conditions ( |3.64j ) to exclude breather-type solutions. 



The sd-CNLS equations can be cast into alternative expressions by some 
transformations. We set a pair of variables and rffl to be constant in (|1.5Q , for 
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instance, 7*1 = 1 and consider a transformation of variables, 

' n c ' n 



Then, we obtain a simplified deformation of (|L 
_ V rj( fc V^ • (a (j) +a {j) 



fc=l 

Q r U) m-1 



j = 1,2, . . . ,m- 1. 



1 ~ / j ' n In V n+1 ^ ' n-ly 

k=l 



To consider a multi-field extension of the full-discrete NLS equation 0, ||, f44j by 
our method, we need to modify the L n -matrix (|2.3|) appropriately. The details of the 
analysis will be reported in a separate paper. 
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Appendix. Trace formulae 

In this appendix, we show interrelations between the generator of the conserved densities 
tr{log(J + Q n T n /z)} for (|2.6|) in section |2]2| and the scattering data in section 3. 



Let us recall that we have transformed the scattering problem (3T) into the 
scattering problem (|3.3j ) by the gauge transformation (|3.2| ). For the moment, we do 
not impose restrictions such as ( |2.31| ) and ( |2.32j ) (or equation ( |3.5| )) on square matrices 
Q n and R n and consider the sd-matrix NLS equation ( |2.6| ). The conditions ( j2.5| ) are 
assumed. We supplement some definitions and relations. First, we define the inverse of 

©by 

Mz) = (f)n{z)A{z) + 4> n {z)B{z) 
1p n (z) = <f> n (z)B(z) - 4> n (z)A{z). 

Secondly, the generator of the conserved densities is invariant under the gauge 
transformation (|3~! 



trlog(/+-g ri $ 2n ^ 
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= trlog(/ + ^Q n ^ 2n ^rn)- (A.l) 
Thirdly, the asymptotic behaviors of the Jost functions <p n and ip n are given by 



(pin 
. ^2n . 




' I ' 




z n 




r\j 


A(z)z n 
_ B(z)z- n _ 


4>ln 
. ^2n _ 




A(z)z n 
_ —B(z)z~ n 



as n — > — oo 



as n — > +00 



as n 



-00 



J 

O 



as n — > +00. 



Further, we can prove that 4>2n4>in * s a polynomial in 1/z and 4>2n4>in is a polynomial 
in z. Therefore, we can replace 4>2n4>in and ^n^fn by r(~) and in section ^]2| 
respectively, except for the difference in gauge (see equation ( |3.2j )). It is important that 
the ratios of two components, 4>2n4>in an d ^in^im are invariant when we consider the 
time-dependent Jost functions <f>$ = n e z2 * and = ip n e z2t . 

Now, we can relate the scattering data with the generator of the conserved densities. 
The determinant of A(z), det A(z), is expressed as 

logdet A(z) = trlogA(z) 



tr [log(0in+i- 



-n— 1 \ 



log(0 ln ^")] 



n=— 00 
00 



tr lo g[01n+101r^ 



1 _-ll 



n=— 00 
00 



1 x 



tr lo g 1 + -Qn(t>2n4>l n 



n=— 00 
00 



tr £ log I+_Q n rH 

L z J 

n=— 00 

°° 1 1 1 
tr X! ["Qn-^-Rn-l-^l H j{QnRn-2 — QnF2R n -lQn-lF2R n -2 ~ -(Qn^Rn-lFl) 2 } + 



(A.2) 



Similarly, det *4(z) is rewritten as 

log det A(z) 

00 

= - tr [log^m+i^"" 1 ) - log(^i^- n ) 



n=— 00 
00 



-tr 51 log[/ + -Qn^nV'ln 
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-tr £ log[/ + -Q n rt +) " 

n=— oo 
oo 

tr [- Mi - QnRn) + 2?Q n F 2 Rn+lF 1 

n=—oo 

+ z {QnRn+2 ~ QnRn+2F\Q n +lRn+lFi 



1 



{QnF2Rn+lF\) 

2 } + 



(A.3) 



Here the time independence of A(z), At{z) = O, is proved in the same manner as in 



section |3.3| . It is now clear how the scattering data are expressed in terms of the integrals 
of motion for the sd-matrix NLS equation. However, it should be stressed that the above 



expansions do not yield local conservation laws. The method presented in section |2T^ is 
useful because it gives not only the densities but also the corresponding fluxes. 

Conversely, the integrals of motion can be expressed in terms of the scattering data. 
For simplicity, we assume that 

i) Q n and R n are expressed as ( |3.41| ). Thus, proposition 3.4. 1 holds. 

ii) det A(z) and det A(z) have AN simple zeros outside and inside the unit circle C, 
respectively. None of them lies on the unit circle C. 

hi) det A(z) and det A(z) approach 1 rapidly as \z\ — > oo and z — > 0, respectively. 
Then, we can derive the following expansion for the sd-CNLS equations Q2.30 ): 

log det A(z) 

AN 



11 1 

E ~n [„ E{ ( 7 * 

n=l Z U j=l Z ] 

°° 1 1 2N 1 

51 Z2k \t 5Z{ (3 

fc=i a ft j=i z j 

00 1 1 ^ 1 
E^fc[^E{(- 



21, 



— z 



w n ~ l log det(A(w)A(w))dw 



2tii Jc 
1 



7T1 JC R 



w ifc_1 log det(A(w) A(w))dw 



2k 



[W 



2k-l 



+ w-^-^hgl det A(w)\ 2 dw . 

(A.4) 



The coefficients of l/z 2k (k = 1, 2, . . .) give an infinite number of the integrals of motion, 
which we call the trace formulae. Here, Cr and Cur denote the right-half portion and 
the upper-right portion of the unit circle C, respectively, as is mentioned in section 
3. It is recalled that z in ( |A.2j ) and z in ( |A.4j ) have the difference 7r/4 in their phases 
(see equation ( |3.4| )). The determinant of A(z) is related to the determinant of A(z) by 
( |3.36| ). It can also be shown that det A(z) and det A(z) are connected by 



log det A(z) = - ^ logdet(7 - Q n R n ) + log det A(z). 

n=— 00 

Thus, we can directly obtain expansions of log det A(z) and logdet^4.(z) with respect to 
z from ( JO! ). 

A derivation of ( |A.4| ) is omitted because it is analogous to that in the continuous 
theory ^ |56], |57| . Related results were also obtained by Kodama [p5| . 
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